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Using density functional calculations we explore the properties of as-yet-unsynthesized 3d — 5d 
ordered double perovskites (A2BB'Oa) with highly polarizable Bi 3+ ions on the A site. We find that 
the Bi2NiReOe and Bi2MnReOe compounds are insulating and exhibit a robust net magnetization 
that persists above room temperature. When the in-plane lattice vectors of the pseudocubic unit 
cell are constrained to be orthogonal (for example, by coherent heteroepitaxy) , the ground states 
are ferroelectric with large polarization and a very large uniaxial magnetocrystalline anisotropy with 
easy axis along the ferroelectric polarization direction. Our results suggest a route to multiferroism 
and electrically controlled magnetization orientation at room temperature. 
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INTRODUCTION 



There is increasing current interest in developing mul- 
tiferroic materials with ferromagnetic and ferroelectric 
order in the same phase for future spintronic or mag- 
netoelectronic devices [3-(3|. Although new multiferroics 
are being predicted and synthesized at an accelerating 
pace, a major obstacle for their adoption in applications 
remains their low magnetic ordering temperatures, which 
are generally far below room temperature. The problem 
more generally lies in the scarcity of insulators with any 
net magnetization - either ferro- or ferrimagnetic - at 
room temperature, even before the additional require- 
ment of polarization is included. Most work on novel 
multiferroics has been restricted to 3d transition metal 
oxides (see, for example, Refs. [BT-flijl). motivated by the 
expectation that 4d or 5d compounds would likely be 
metallic. However . recent work on new 3d — 5c? double 
perovskites [I4j showed this to be a misconception, 
and identified a ferrimagnetic insulator, S^CrOsOg, with 
magnetic ordering temperature well above room tem- 
perature. First-principles calculations have been shown 
to accurately reproduce the measured magnetic ordering 
temperatures, 15, 1(| in the Sr 2 CrM06 series (M = W, 
Re, Os), [13] and have been invaluable in explaining the 
origin of the robust ferrimagnetic ordering. (l8l - [2l| 

We build here on these recent developments to pro- 
pose a route to achieving multiferroics with higher mag- 
netic ordering temperatures. We use 3d — bd dou- 
ble perovskites to achieve high temperature magnetism 
combined with insulating behavior, and introduce lone- 
pair active cations on the A sites to induce polarizabil- 
ity [22I, ■ The heavy Bi and 5d elements provide an 
additional desirable feature: The spin-orbit interaction 
strongly couples the magnetic easy axis to the ferroelec- 
tric polarization direction. Specifically, we explore two 
compounds, Bi 2 MnRe06 and Bi 2 NiRe06- We use den- 
sity functional theory to calculate their zero Kelvin struc- 



ture and magnetic ordering, and Monte Carlo simulations 
with parameters extracted from the first-principles cal- 
culations to calculate their magnetic ordering tempera- 
tures. We find that both materials are magnetic insula- 
tors, with ordering temperatures well above room tem- 
perature. While the global ground state in both materi- 
als is anti-polar, we show that when the in-plane lat- 
tice vectors form an angle close to 90° a ferroelectric 
phase results. In practice this could be achieved using 
coherent heteroepitaxy and the resulting ferroelectric- 
ity could be manipulated using strain. We show that 
the magnetic behavior results from an antiferromagnetic 
superexchange ujj, [2J] and the ferroelectricity from the 
Bi 3+ lone pairs, reminiscent of other Bi-based multifer- 
roics. Finally, we investigate the effect of strong spin- 
orbit coupling (SOC) due to the presence of heavy Bi 3+ 
and Re 4+ cations on the structural, magnetic, and ferro- 
electric properties of these compounds. 



METHOD OF CALCULATION 

Structural optimizations were performed using the Vi- 
enna Ab-initio Simulation Package (VASP) with projec- 
tor augmented wave (PAW) potentials; [25j the scmicore 
p states of Mn, Ni, and Re and the d states of Bi were 
included in the valence. We used an energy cutoff of 500 
eV and 6x6x6, 4x4x6, and 4x4x4 T-point centered 
fc-point meshes for unit cells containing 10, 20, and 40 
atoms respectively. The exchange-correlation functional 
was treated within the generalized gradient approxima- 
tion (GGA); [H] while extension to the GGA+U method 
had only a small influence on the electronic properties, we 
discuss later its effect on the structural behavior. Mag- 
netic ordering temperatures were obtained using a finite- 
temperature Monte Carlo scheme within a Heisenberg 
Hamiltonian H = — ^ V\ - Ji ,-rVL • M, , where M< and 



Mj are the magnetic moments on sites i and j of the 
crystal lattice. The exchange constants Ji.j were cal- 
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culated in the frozen-magnon scheme 27J using the all- 
electron full-potential linearized augmented plane-wave 
code FLEUR (28|, with a plane- wave cutoff of 4.2 hartrees 
and 15x15x15 /c-points and a 6 x 6 x 6 spin-spiral 
grid in a 10-atom unit cell. Ferroelectric polarizations 
were extracted from the shifts of the centers of Wannier 
functions (29j . The magnetocrystalline anisotropy energy 
was extracted from the self-consistent total-energy calcu- 
lations within the FLEUR code, on a 13 x 13 x 13 fc-point 
grid. 

First, we calculate the lowest energy structure for both 
Bi 2 MnRe06 and Bi 2 NiReOg within the GGA approxi- 
mation. We proceed by successively freezing in the 12 
centrosymmetric combinations of tilts and rotations of 
the oxygen octahedra that can occur in ordered double 
perovskites [30] . We then further reduce the symmetry of 
these lowest energy tilt systems by displacing the anions 
and cations relative to each other in the manner of a po- 
lar displacement. Next we optimize the atomic positions 
(by minimizing the Hellmann-Feynman forces), as well as 
the unit cell shape and volume within each symmetry. In 
this first series of optimizations we assume ferrimagnetic 
ordering and do not include SOC; the influence of SOC 
on the lowest-energy structures is examined later. 



GROUND-STATE STRUCTURE AND 
INFLUENCE OF EPITAXIAL CONSTRAINTS 

Our calculations yield the same monoclinic, cen- 
trosymmetric P2\/n symmetry ground state for both 
compounds. It is characterized by an a~a~c + tilt pat- 
tern of the oxygen octahedra. Two lattice vectors are 
identical with angles of 93.6° and 93.5° between them 
in the Mn and Ni compounds, respectively (in the fol- 
lowing we refer to these as the "in-plane" lattice vec- 
tors); the third lattice vector is perpendicular to the first 
two and of different length. Importantly, in both com- 
pounds there is a slightly higher energy ferroelectric (FE) 
structure with rhombohedral i?3 symmetry (see Table [TJ) . 
The i?3 structure corresponds to an a~a~a~ tilt pattern 
yielding lattice vectors of equal length and rhombohedral 
angles of 60° and 61° for the Mn and Ni compounds re- 
spectively, combined with polar relative displacements of 
anions and cations along the [111] rhombohedral axis in- 
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Bi 2 MnRe0 6 


Bi 2 NiReO 
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TABLE I: Total energies (per formula unit) of the lowest- 
energy phases of Bi2MnRe06 and Bi 2 NiRe06 and the cubic 
Fm3m phase with respect to their ground state P2\/n phase. 
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FIG. 1: (Color online) (a) Ground-state P2\/n structure and 
(b) the ferroelectric R3 structure of the double perovskites 
Bi 2 MnRe0 6 and Bi 2 NiReOe. Bi 3+ cations are depicted as 
the large (brown) spheres and O 2- as the small (red) spheres. 
Re 4+ and Mn 2+ (Ni 2+ ) are alternating in the shaded octa- 
hedra, in a three-dimensional checkerboard manner (this is 
indicated by the color of the octahedra) . While in the P2\ jn 
symmetry the Bi 3+ cations form an "antipolar" pattern, their 
coherent shift along the [111] direction can be clearly seen in 
the R3 structure. 



duced by the well-established stereochemically active Bi 
lone pair [22| (see Fig. [1]). The energy of the FE phase 
is significantly lower than that of the corresponding cen- 
trosymmetric i?3 phase, suggesting a high ferroelectric 
ordering temperature if this symmetry could be stabi- 
lized. For comparison, in the prototypical multiferroic 
BiFeOa with Tq of 1103 K, we obtain a difference of 515 
meV per two formula units between the corresponding 
states. The calculated polarization in both cases is along 
the [111] direction, and comparable in size to that of 
BiFe0 3 (see Table [TTJ) . 

While the ground state P2i/n phase is of interest in 
its own right as a possible high temperature magnetic 
insulator, we next explore whether it is possible to iden- 
tify conditions that stabilize the i?3 ferroelectric phase. 
We use the fact that in the FE i?3 phase all three lattice 
vectors have equal length, and the rhombohedral angle is 
equal or close to the ideal value of 60°. In contrast, in the 
P2i/n phase, only the lengths of the two in-plane lattice 
vectors are equal, and the inter-in-plane angles deviate 
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FIG. 2: (Color online) (a) Distorted Bi2MnReOs and Bi2NiRe06 obtained from the R3 (red circles) and P2i/n (blue triangles) 
structure by constraining the two in-plane lattice vectors to be equal and the angle between them to 90°. The out-of-plane 
parameter and the ionic positions were relaxed. The energies are positioned with respect to the bulk unconstrained P2\jn 
ground state (blue squares). The red diamonds show the energies for the unconstrained R3 structure. A crossover between the 
antipolar (AP) and the ferroelectric (FE) state occurs at 1.8% (with respect to the minimum of the FE curve) of compressive 
strain in Bi2MnRe06- Note that the energy of the R3 state lies on the strain curve for Bi2MnReC>6, but is somewhat lower for 
Bi2NiReOs, due to the choice of a square in-plane lattice: the rhombohedral angle in R3 phase of Bi2MnReOe is 60°, while in 
Bi2NiReC>6 it is 61°. (b) Similar to (a), with an addition of a Hubbard U on Mn (3 eV) and Ni (6 eV) 3d states. The shaded 
area shows the AP-FE crossover region. 



strongly from the ideal 90°. Motivated by these observa- 
tions, we first constrain all three lattice parameters to be 
equal in length and indeed find that this constraint sta- 
bilizes the FE phase. Next we enforce that only two of 
the three lattice parameters are equal in length, but con- 
strain them to be perpendicular to each other and again 
find that the FE phase is the lowest energy 37j . Such a 
constraint could be achieved in practice through coherent 
heteroepitaxy between a thin film and a substrate, and 
is often described as a biaxial strain state in the litera- 
ture. In Fig. Uta), we show the calculated energies of the 
previously i?3 and P2\/n structure types subject to this 
additional constraint, with the in-plane lattice parame- 
ters varied over a range of realistic substrate strains corre- 
sponding to (001) epitaxial growth. For each in-plane lat- 
tice parameter we allow the length and angle of the out- 
of-plane lattice parameter and atomic positions to relax 
to their lowest-energy configurations. Our main finding 
is that, under the investigated constraint, the FE phase 
is lower in energy than the antipolar (AP) phase. At the 
optimized value of the in-plane lattice parameter (8.16 A 
for Bi 2 MnRe0 6 and 8.07 A for Bi 2 NiRe0 6 ) the energy 
differences between the FE and AP phases are 36 and 
30 meV per formula unit, respectively. The FE phases 
are robustly insulating with calculated band gaps of 0.7 
and 0.3 eV, respectively, and have large polarizations of 
84 /xC/cm 2 (Bi 2 MnRe0 6 ) and 81 ^C/cm 2 (Bi 2 NiRe0 6 ) 
along the pseudocubic [111] direction. Notice, however, 
that as in-plane compressive strain is applied, the dif- 
ferences in energies reduce, and in fact in Bi 2 MnRe06 a 
crossover is reached at a compressive strain value of 1.8%. 



We therefore expect that both compounds, although AP 
in bulk, will in fact be ferroelectric and hence multiferroic 
in thin-film form, over a range of experimentally accessi- 
ble strains. We see later that SOC tips the scale toward 
the FE phase, while electron correlations favor the AP 
phase. 



ELECTRONIC AND MAGNETIC PROPERTIES 
OF THE FERROELECTRIC PHASE 



Next, we analyze the magnetic behavior of the two 
ferroelectric compounds. We find that in both cases 
the lowest-energy ordering is ferrimagnetic, with the mo- 
ments on the 3c? and 5c? transition metal ions antialigned 
(see Fig. [3]). In both compounds, Re is in formal oxida- 
tion state 4 + corresponding to a filled d — t 2g manifold in 
the minority spin-channel. The oxidation states of Mn in 
Bi 2 MnRe0 6 and of Ni in Bi 2 NiRe0 6 are 2+: t\ g e 2 (high 
spin) for Mn and i 2g e 2 for Ni. Due to these specific con- 
figurations of the outer electronic shells, and the nearly 
150° Mn-O-Re (Ni-O-Re) angle, the magnetic moments 
on Re and Mn/Ni are coupled via an antiferromagnetic 
superexchange mechanism. {33T435T ] The result is a ferri- 
magnetic configuration with a total spin moment of 2 [Ib 
in Bi 2 MnReC>6 and 1 /is in Bi 2 NiReOg. 

Our calculated magnetic ordering temperatures in the 
ferroelectric phase are 330 K for Bi 2 MnRe06 and 360 K 
for Bi 2 NiRe06, both significantly above room temper- 
ature. As a comparison, using the same method, we 
calculate a magnetic ordering temperature of 670 K for 
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FIG. 3: (Color online) Density of states (DOS) of the ferro- 
electric phase of Bi 2 MnRe0 6 (a) and Bi 2 NiReO e (b). The 
local magnetic moments on Mn/Ni and Re are antialigned. 
Thin black line, total DOS; green dashed line, Re; orange 
solid line, Mn/Ni. 



S^CrOsOe (the experimental value is 725 K). Note that 
if growth can only be achieved in the form of ultra-thin 
films, the ordering temperatures will likely be reduced. 



SPIN-ORBIT COUPLING AND CORRELATIONS 
EFFECTS 

Since the heavy Re and Bi atoms in our double per- 
ovskites are also likely to exhibit strong SOC, and con- 
sequently significant magnetostructural coupling, |18[ we 
repeat our calculations with SOC explicitly included. We 
find that the ground state remains P2\/n AP, but now its 
energy difference from the FE state with i?3 symmetry 
is reduced to 10 meV in P^MnReOg and only 2 meV in 
Bi2NiReOg. Both compounds still have semiconducting 
gaps which are now somewhat reduced while the ferro- 
electric polarization of the i?3 phase is slightly increased 
(see Table HI)) . The total magnetic moments (spin + or- 
bital) amount to 2.34 and 0.58 hb, respectively. Note 
that, in contrast to half-metallic ferromagnets where the 
SOC introduces states in the minority-spin gap yielding 
a noninteger total magnetic moment, here the noninte- 
ger moment is the consequence of the mixing of spin-up 
and spin-down states (spin is no longer a good quantum 
number with SOC included), while the gap is preserved. 

In both compounds the magnetic easy axis lies along 
the ferroelectric polarization direction. The associated 
anisotropy energy is large and comparable to that in cur- 



TABLE II: Properties of the bulk ferroelectric R3 phase in 
Bi 2 MnRe0 6 and Bi 2 NiRe0 6 with and without the SOC in- 
cluded: band gap, ferroelectric polarization P and the total 
magnetic moment M (spin+orbital in case of SOC) per for- 
mula unit (f.u.). 



rent magnetic recording media: it amounts to 7 and 5.5 
meV per formula unit in Bi2MnReOg and Bi2NiRe06, 
respectively. This suggests an exciting possibility of elec- 
trical control of the magnetization direction. 



Finally, we investigate the influence of electronic cor- 
relations by adding a Hubbard U on the d states of Mn 
and Ni. In order to reduce the computational effort, we 
do not include the SOC. We keep in mind, however, that 
a consequence of its inclusion, as we have seen, is a re- 
duction in the energy of the FE state with respect to 
the AP one. Introduction of Hubbard U corrections (on 
Mn U=3 cV and J=0.87 eV, and on Ni U=6 eV and 
J=0.90 eV) does not change the ground state structure 
{P2\/n still has the lowest energy). Interesting quanti- 
tative changes occur in the strain dependence however 
[Fig. E{b)]. While within the GGA the ferroelectricity 
was robust to the choice of substrate lattice parameters 
over a likely range of accessible strains provided that the 
in-plane lattice parameters were constrained to be per- 
pendicular to each other, in GGA+C7 we find a cross-over 
between the FE and AP states at moderate strain values. 
This suggests the intriguing possibility that an external 
strain, for example from a piezoelectric substrate [36[ 
could induce a "dipole-flop" transition from a non-polar 
structure in which the dipoles are antialigned to a ferro- 
electric one. For the assumed value of U, the crossover 
region [shaded area in Figdfb)] in Bi2MnReOg is within 
less than 1% from either of the two minima: a small 
compressive strain will push the system toward the AP 
state, while a small tensile strain will favor a FE state. 
In Bi2NiRe06 the applied U pushes the AP state lower 
in energy, but the crossover region to the FE state is still 
within an experimentally accessible range of 1.5% tensile 
strain. While we can not make a quantitative prediction, 
since the relative energies of the AP and the FE state de- 
pend strongly on the value of the applied U correction, 
within a reasonable range of U the crossover region lies 
within 2% strain from the ground state structure. 
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CONCLUSION 

In summary, we predict from first-principles calcu- 
lations that the 3d — 5d double-perovskite compounds 
Bi2MnRe06 and P^NiReOg are high-Tc ferrimagnetic 
insulators. Moreover, it is likely possible to stabilize them 
in thin-film form in a ferroelectric phase with high polar- 
ization along the [111] direction. The magnetic easy axis 
lies along the ferroelectric polarization direction with a 
very large associated anisotropy energy. The different 
sizes and charges of the 3d and 5d transition metal ions 
suggest that the required -B-site ordering should be ex- 
perimentally feasible; we hope that our findings will ini- 
tiate experimental efforts to synthesize these novel mul- 
tiferroic compounds. 
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